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The full spectrum of two interacting electrons in a disor- 
dered mesoscopic one-dimensional ring threaded by a mag- 
netic flux is calculated numerically. For ring sizes far ex- 
ceeding the one-particle localization length L\ we find sev- 
eral /t/2e-periodic states whose eigenfunctions exhibit a pair- 
ing effect. This represents the first direct observation of 
interaction-assisted coherent pair propagation, the pair be- 
ing delocalized on the scale of the whole ring. 

PACS numbers: 72.15, 73.20 

For more than three decades Anderson localization has 
been a subject of intensive research (for a review see [1]). 
Among the first results in this field was the fact that non- 
interacting electrons in one dimension (Id) are always 
localized and that their wave-functions decay exponen- 
tially on the scale of the localization length L\. Very 
recently, Shepelyansky [2] and later Imry [3] considered 
two interacting electrons in a Id random potential. Both 
authors suggest that in a regime of not too strong dis- 
order correlated electron pair states may extend over a 
distance L'j oc L\ much larger than L\. This completely 
novel effect could very well have a fundamental impact on 
our understanding of the localization problem and other 
related areas. However, the methods employed in [2,3], 
though powerful and suggestive, are partly approximate 
and partly qualitative. Shepelyansky's original reason- 
ing relied on certain assumptions of a statistical nature, 
allowing him to map the problem onto a random band 
matrix model. A numerical study of this model together 
with additional evidence from a variety of other models 
formed the basis for his claim. Imry, on the other hand, 
invoked the Thouless scaling block picture to arrive at 
precisely the same results as Shepelyansky. Moreover, 
Imry's approach seems to be very suitable to generalize 
the effect to higher dimensions. 

The key quantity in [2,3] was the distribution of in- 
teraction matrix elements in the disorder-diagonal basis. 
This distribution defines the statistical properties of the 
band matrix model and provides the interblock coupling 
in the scaling approach. In a recent study [4] employing 
the transfer matrix technique we investigated numerically 
both finite and infinite systems starting from first prin- 
ciples. In this work, the principal effect was confirmed, 
but the two-particle localization length was found to 
scale with a smaller exponent, L'j oc L\ 65 , in the param- 
eter range investigated. This modification could be at- 
tributed to particular properties of the above-mentioned 



distribution which had not been taken into account in 
[2,3]. 

In the present paper, we demonstrate the existence 
of /i/2e-periodic states in the spectrum of a disordered 
mesoscopic ring with two interacting electrons by direct 
diagonalization. Inspection of the two-electron eigen- 
functions reveals that this period halving is indeed due 
to a pairing effect: Both electrons propagate coher- 
ently around the ring, staying within a distance of a 
few L\ from each other. This is the first direct ob- 
servation of interaction-assisted coherent pair propaga- 
tion in disordered systems. The fact that the effect can 
be observed for the ring sizes considered (circumference 
N = 50,100,150) and the disorder chosen (Li k, 11) 
supports our earlier investigations [4], where we found 
L2 ~ 25 > L\ for these system parameters and U = 1 
(see (1) below). Here and in the remainder of the text we 
set the lattice spacing to unity. The flux sensitivities of 
the h/2e states agree with what one can expect: an expo- 
nential reduction exp(— N/L2) for N > L2. Having cal- 
culated the full spectrum of the two-electron ring we also 
comment on the additional observation of interaction- 
enhanced /i/e-oscillations. In the less localized case the 
enhancement simply originates from a screening effect, 
while sensitive h/e states in the large rings benefit from a 
partial pairing effect. Our findings could also be relevant 
for the persistent current problem (see [5] and references 
therein). 

To be specific, we consider a Id ring which is threaded 
by the magnetic flux <p = $/$0; where $0 = h/e is the 
flux quantum. We employ a tight binding model with 
N sites and random on-site energies V n (n = 1, . . . , N), 
the latter being uniformly distributed in the interval 
[—VP, VP]. The electron-electron interaction is described 
by a local Hubbard term with strength parameter U . 
Thus the Hamiltonian reads 
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The operators c+ a and c n a create and destroy an elec- 
tron at site n (we set c^j^ a = c^) with spin a =|, J,, 
respectively. 

We restrict our investigation to the case of two elec- 
trons with opposite spins. Then, the spatial part of 
the wave-function is symmetric, double-occupancy is al- 
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lowed and the Hubbard interaction is relevant. We work 
in the basis of the M = N(N + l)/2 different wave- 
functions 

I ^ lt n 2 ) \cr) = -j= (c+ T c+ d + c+ ^c+J I 0} . (2) 

Since the Hamiltonian does not couple the spin compo- 
nents, we drop the antisymmetric spin part | a) of the 
wave-function. 

The resulting M x M Hamiltonian matrix is diagonal- 
ized by means of the Lanczos algorithm. The sensitivity 
of the two-electron energies E m (m = 1, . . . , M) on the 
flux as well as their periodicity is determined by calcu- 
lating the full spectrum of M eigenvalues of % for a few 
values of the magnetic flux between <p = and <p = 1/2. 
For some levels with interesting flux dependence we also 
determine the corresponding eigenfunctions |\? m ). The 
shape of these eigenfunctions, i.e. the spatial distribu- 
tion of the coefficients ^m{n\,n 2 ) = (ipn lt n 2 I ^m), is 
very useful for the interpretation of the properties of the 
states and also helps to clarify the effects which lie at the 
origin of their behavior. Furthermore, we determine the 
"local" current 

I m (ni,n 2 ) = + «2)* m («i, n 2 ) (3) 

- *m(«i, «2)* m («i + 1, n 2 )} + {ni <-> n 2 }) , 

corresponding to the center-of-mass motion of the two 
electrons. This quantity yields additional information 
about the way in which the two electrons propagate. The 
total current carried by the state | \? m ) is then given by 

We briefly discuss the system parameters used in the 
sequel. Present-day computer technology restricts us to 
ring sizes of the order of N ps 10 2 . We have calculated full 
spectra for N = 50, 100, and 150. The disorder parame- 
ter W has to be chosen such that the ring is strongly lo- 
calized in order to unambiguously identify coherent pair 
propagation. On the other hand, to see an effect at all 
we must have a sizable enhancement factor L 2 /L\ for 
the pair localization length L 2 . We choose W = 1.5 and 
hence L\ = 25/ W 2 ps 11 in the present work. For this 
value (and U = 1) we found L 2 ps 25 in our transfer ma- 
trix study [4], i.e. an enhancement factor of more than 
two. For the Hubbard interaction we consider the values 
U = —1, 0, 1, as in [4]. We can estimate the amplitude A 
of the h/2e oscillations by the following considerations. 
For a ring of size L 2 roughly LiL 2 states will be cou- 
pled by the interaction, since a pair state has width L\ 
and length L 2 . The amplitude A at scale L 2 is given by 
some fraction c of the level spacing A 2 (L 2 ) of the cou- 
pled states. The total bandwidth being 2(4+ W), one gets 
A ps 2c(4+W)/LiL 2 . At scales N > L 2 this amplitude is 
suppressed by exp(— N/L 2 ). With W = 1.5 we therefore 
expect A ps c/200, c/1500, c/10000 for the most sensitive 
states in rings with N = 50, 100, 150, respectively. 



In Fig. 1 we compare the flux sensitivities of the en- 
ergy levels for different ring sizes and different values of 
U . We have plotted the difference between the maximum 
and the minimum value in the interval < <p < 1/2 (i.e. 
the amplitude A) versus the energy value at zero flux. 
As the ring size increases the overall amplitudes decrease 
exponentially. More and more levels become localized 
and thus independent of the flux. Investigating the pe- 
riodicity of the flux-sensitive levels, one finds for U ^ 
that levels with strong h/2e admixtures play a more and 
more prominent role in larger rings. These admixtures 
can be identified by the non-monotonous behavior of the 
eigenvalue in the flux interval < <p < 1/2 [6]. The am- 
plitudes of the most sensitive h/2e states are indicated 
in Fig. 1 by horizontal lines. For N = 50 the separation 
of the electrons forming a pair can be of the order of the 
ring circumference so that one expects h/e rather than 
h/2e oscillations. Indeed, very few h/2e states exist for 
N = 50 and their amplitudes are dominated by those 
of the h/e states by orders of magnitude. Most strik- 
ingly, this is no longer true for N = 100. There, h/2e- 
oscillations are less than an order of magnitude smaller 
than those of the largest /i/e-periodic states. With an 
amplitude of about 10 -4 the h/2e states range between 
the comparatively few states with appreciable flux sensi- 
tivity, demonstrating the growing importance of the pair 
propagation process. The trend towards ever more pro- 
nounced /i/2e-oscillations continues for N = 150. There, 
only very few of the M = 11325 levels are more sensi- 
tive to the flux than the /i/2e-periodic ones. Assuming 
that some of the h/e states for N = 50 are actually pair 
states we can compare our estimates for the amplitude 
A with the numerical results. The exponential factor 
exTp(—N/L 2 ) (and not exp(— N / L\)) is clearly confirmed 
by the data. Moreover, the absolute values are in good 
agreement for c ps 1/6. 

In Fig. 2 we show the flux-dependence for two se- 
lected h/2e states, one for N = 100 and the other for 
N = 150. Density plots of the amplitudes |\? m (ni, n 2 )\ 2 
of the corresponding wave-functions and the local cur- 
rents I m (ni, n 2 ) are presented in Fig. 3 and Fig. 4, respec- 
tively. They clearly exhibit the "cigar-shape" character- 
istic for coherent pair propagation: The wave-functions 
are concentrated around the diagonal n\ = n 2 with a 
transverse extension given by a few L\, but in the lon- 
gitudinal (diagonal) direction they extend over the full 
ring. We recall that the circumference of the two rings 
is given by 9_Li and 1AL\, respectively. The strong h/2e 
admixture to the flux dependence of the eigenvalues to- 
gether with the shape of the wave-function convincingly 
demonstrates that two electrons, staying within a dis- 
tance L\ of each other, propagate as a composite entity 
with charge 2e. Moreover, the local current confirms and 
reinforces this interpretation. For instance, the branch 
along n ps 60 in Fig. 4 is clearly identified as a dead 
end, showing that the electrons do not propagate around 
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the ring independently. This is consistent with the al- 
most perfect /i/2e-periodic behavior of the correspond- 
ing energy level (see Fig. 2). We emphasize that the local 
current defined in (3) measures only the center-of-mass 
motion of the electrons. This quantity was chosen to 
simplify the presentation. We have checked that the local 
currents associated with the relative motion are similarly 
concentrated around the diagonal. 

The behavior of the /i/e-periodic states is also of inter- 
est. Comparing the flux-sensitivities in Fig. 1 for U = 
and U = 1 one notices that while the average amplitude 
is slightly reduced for U ^ 0, a few very sensitive states 
emerge. This has been found for repulsive (U = 1) as well 
as for attractive (U = — 1) interactions. This interaction- 
enhanced electron mobility is of a less subtle nature than 
the pair propagation effect. In Fig. 5 (left) we show a 
density plot of the wave-functions of particularly flux- 
sensitive /i/e-periodic states for U = 1 and N = 50. 
This state exhibits a very pronounced cross-like struc- 
ture, proving that the electrons propagate independently. 
The enhancement is simply due to screening: One of the 
electrons fills the deepest hole in the random potential 
landscape, occupying a state similar to the one-particle 
ground state |xo}- The second electron then moves in the 
smoothed effective potential = V n + |xn| 2 , where \ n 

is the component of |xo) associated with site n. There- 
fore, this electron becomes delocalized and very sensitive 
to the magnetic flux. Similarly, for U = — 1 the strongest 
positive fluctuation of the random potential is reduced 
by the attractive interaction between the delocalized elec- 
tron and its localized partner, the latter occupying a one- 
electron state of high energy. This explains why highly 
flux-sensitive h/e states of this type tend to occur in the 
upper (lower) half of the spectrum for U = —1 (U = 1). 
On the right hand side of Fig. 5, we show a particularly 
flux sensitive h/e state for N = 100 and U = 1. Here, 
the comparatively high amplitude benefits from the pair- 
ing effect discussed above in the context of /i/2e-periodic 
states. For this state, the pair is not extended over the 
whole ring, but splits up before each electron completes a 
cycle. The contribution of pair propagation in h/e states 
for N = 150 is even more pronounced. 

In conclusion, we have found /i/2e-periodic states in 
the spectrum of a Id disordered mesoscopic ring with 
two interacting electrons and circumference N > L\. 
The corresponding wave-functions and local currents un- 
ambiguously demonstrate the pairing effect proposed by 
Shepelyansky. The amplitudes of these states agree well 
with simple estimates. The method, direct diagonaliza- 
tion, is free from any approximations whatsoever. The 
present investigation is in good agreement with our pre- 
vious quantitative study of the pair localization length 
L2 using the transfer matrix [4]. Furthermore we have 
observed an interaction-induced enhancement of h/e— 
oscillations due to the screening of the disorder potential 
in small rings, while pair propagation effects play an im- 



portant role even for the /i/e-oscillations in larger rings. 

One might argue that the pair propagation effect is of 
limited relevance since we found only a few flux-sensitive 
states somewhere in the spectrum of the mesoscopic ring. 
This is, however, due to the relatively small system sizes 
to which we are restricted. For N ps L\ 1 

the effects described in this paper should be much more 
dramatic. Furthermore, in a regime where the one- 
particle localization length L\ at the edges of the one- 
particle band becomes much larger than the lattice spac- 
ing one can expect even the lowest-lying states of the 
ring to be affected. Also, both coherent pair propagation 
and the mechanism underlying the enhancement of h/e— 
oscillations might be relevant for the persistent current 
problem. In this case, one deals with a diffusive ring so 
that the maximum pair separation is larger than the ring 
size and one should not expect /i/2e-oscillations. It is 
quite attractive to speculate that coherent pair propaga- 
tion and the effective weakening of the disorder potential 
proposed in [5] are manifestations of the same physical 
effect at different scales. Much further work is needed to 
clarify these issues and to extend the results obtained up 
to now to higher dimensions and higher electron number. 
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FIG. 1. The amplitudes A of the oscillations of all levels 
as a function of the energy of the states. The upper pictures 
show the cases U = (left) and U = 1 (right) for N = 50. 
The lower ones correspond to N = 100 (left) and N = 150 
(right) at U = 1 (note the different scales for A). The total 
number of levels M is 1275, 5050, and 11325 for the ring sizes 
50, 100, and 150, respectively. The horizontal lines indicate 
the most sensitive /t/2e-periodic levels. 
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FIG. 3. The state with the highest h/2e flux sensitivity at 
N = 100 and U = 1 corresponding to the level shown in Fig. 2 
(top) for the flux value tp = 0.15. The absolute square of the 
wave-function is shown on the left-hand side. Dark regions 
correspond to high values. The associated local current is 
also shown (right). Here, dark regions correspond to positive 
values (propagation towards increasing n\ and/or n^), while 
bright regions indicate the opposite sign. In the grey areas, 
the current is close to zero. 
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FIG. 4. A state with a high h/2e flux sensitivity at N = 150 
and U = 1 corresponding to the level shown in Fig. 2 (bot- 
tom). The absolute square of the wave-function (left) and 
the local current (right) are shown for tp = 0.15. 
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FIG. 2. The flux dependence of /t/2e-periodic levels with 
high sensitivity for a ring of size N = 100 (top) and N = 150 
(bottom) at U = 1. While the first one still has a small 
admixture of an /t/e-periodic oscillation, the curve for the 
larger system exhibits almost perfect /t/2e-periodicity. 
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FIG. 5. The states with the highest flux sensitivities in 
rings containing N = 50 (left) and N = 100 (right) sites. 
They are found at E ~ —2.411 and E K, —1.139, respectively. 
In the latter case a tendency towards a concentration near the 
diagonal indicates the formation of a pair, but the electrons 
split up before they complete a cycle. 
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